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Abstract 

For quantum (quasi)particles living on complex toboggan-shaped curves which spread 
over N Riemann sheets the approximate evaluation of topology-controlled bound- 
state energies is shown feasible. In a cubic-oscillator model the low-lying spectrum 
is shown decreasing with winding number N . 
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1 Introduction 



In one of the most up-to-date reviews of the so called VT — symmetric quantum 
mechanics [1] the existence has been emphasized of multiple connections between 
certain less usual differential Schrodinger equations for bound states and a number 
of important integrable models in statistical physics and/or in conformal quantum 
field theory [2] . This correspondence offers one of explanations of the current growth 
of interest in many apparently exotic bound-state problems 

(f n (q) = E n (p n (q) (1) 

where one considers, say [3], the three-parametric family of complex interactions 

V eff (g) = g 2 (zg) 2M " 2 - a + , M > 1 . (2) 

Apparently, this manifestly non-Hermitian model can hardly find any direct applica- 
bility in quantum phenomenology Fortunately, this first impression proved wrong. 
The relevance of similar non-self adjoint models has recently been revealed not only 
in quantum mechanics (recall, e.g., the so called interacting boson models in nuclear 
physics but also in certain quantum field theories [5] and even in non-quantum 
nonlinear optics |6] (cf. also the other recent applications of non-Hermitian Hamil- 
tonians collected in review papers [7), [8J, [9] and in proceedings [TO]). 

Even though the study of similar models is quite common in mathematics [TT] . 
virtually all of the above-mentioned dedicated reviews of their role in physics only 
appeared very recently. Indeed, the manifestly non-selfadjoint problems of the form 
of Eq. ([I]) + (0) were usually considered "unphysical" . Their studies only appeared 
in the purely methodical context of perturbation theory where, typically (cf., e.g. 
[T2]). differential Eq. ([1]) + ([2]) with M < 2 has been assumed integrated along the 
straight real line or along its shifted, complexified version 

q = s — ie = q^(s) , s G (— oo, oo) . (3) 

A decisive progress has been achieved when, in Ref. [3] , the reality of the spectra 
generated by potentials ([2]) has been given a rather nontrivial proof requiring just 
the following very elementary sufficient condition 

M+ 1 + \2£ + 1| > a. (4) 

One appreciates, in particular, that for any coupling a and dominant exponent 2M 
this reality of the energies (i.e., in principle, their observability) will always be guar- 
anteed once we select a sufficiently large strength £ ^> 1 of the singularity. 

The latter observation attracted our attention for two reasons. Firstly, we kept 
in mind that the presence of the singularity in the potential may imply that the 
individual wave functions (p n (<l) become multisheeted when treated as analytically 
continued functions of complex "coordinate" q [T3]. In parallel, we imagined that the 
condition of physical acceptability (jlj) coincides with the condition of mathematical 
consistency of the so called l/£ perturbation expansions (for details and further 
references see Ref. [H] or Appendix A below). 
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In what follows we intend to make full use of the latter unexpected coincidence 
of the physical and mathematical appeal of the "strongly spiked", t ^> 1 versions 
of potentials (121). An immediate motivation of our present return to the related 
Schrodinger Eq. (JTJ arose from Refs. [T31 [To]. There, the topologically trivial in- 
tegration path fl3]) has been replaced by the complex curves which were allowed to 
encircle the singularities of the potentials (cf., e.g., Refs. [16] and [17] or Sec. [2]below 
for a compact introduction to this new possibility of quantum model-building). 

Unfortunately, even in the simplest models ([I]) the tobogganic bound-state prob- 
lems were found extraordinarily hard to solve numerically [HJ [19]. For this reason, 
no sufficiently reliable illustration of a tobogganic spectrum is at our disposal at 
present. At the same time, the above-mentioned tests of applicability of \ ji pertur- 
bation expansions to non-tobogganic non-Hermitian models [H] may be perceived 
as a source of new optimism and as a decisive encouragement of our present project 
of an approximative evaluation of the tobogganic low-lying spectra. 

Our results will solely involve tobogganic versions of the concrete toy model ([2]) 
possessing the single singularity in the origin q = and considered just at a fixed 
sample exponent M = 3/2 and at the simplest coupling a = 0. Under these assump- 
tions our main attention will be paid to our Schrodinger Eq. ([1]) integrated along 
the whole family of nontrivial, "tobogganic" curves denoted by the symbol q^ N \s) 
and classified, in an exhaustive manner, by their winding number iV = 0, 1, 2, . . . (cf. 
their sample in Figured]). By definition (cf. Sections [2] and [3] below) these smooth 
complex curves will interconnect an N— plet of sheets of the Riemann surface 1Z 
supporting the analytic, multivalued wave functions <p n {q)- 
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Figure 1: The complex straight contour q^°\s) of eq. ([3]) (thin line) and its first two 
winding descendants q^(s) (medium curve) and q^(s) (thick curve). 

We shall reflect the dominance of the centrifugal repulsion near the origin as an 
encouragement of the use of the approximative large— £ method of the solution of 
Eq. ([1]). The essence of this method will be explained in Section HI For our M = 3/2 
(i.e., imaginary cubic) quantum toboggan the application of this method will finally 
be described in Section [5] complemented by a few comments and a brief summary in 
Section [6j 



2 Topologically nontrivial integration paths 



Several important theoretical reasons of the contemporary growing interest in non- 
tobogganic as well as tobogganic models sampled by Eqs. (OQ) + (j2J) + (EJ) deserve 
to be re-emphasized. One of the most typical ones lies in the emergence of several 
unexpected conceptual links between the ODE language and its correspondence to 
the physical theories of integrable statistical and field-theoretical models. In this 
sense, in a way mentioned in (TJ [3], the integrable SUSY point of the sine-Gordon 
model has been found related to the value of the exponent M = 1/2 in (El). Similarly, 
the partnership with the YangLee model emerged at M — 3/2, etc (cf. section 5.6 
in [7J for a "full dictionary"). 

One should also mention the importance of the flow of inspiration in the opposite 
direction. Indeed, certain very specific questions asked in the context of conformal 
field theory [2] brought, as their direct consequence, the immediate turn of attention 
from the most elementary and regular £ = versions of Schrodinger Eq. (j2J) to 
its specific singular generalizations. A concrete example of the emerging family of 
relevant generic potentials with multiple centrifugal-type singularities may be found, 
e.g., as Eq. (1) in Ref. [2]. 

Naturally, after similar generalizations and incorporation of singularities even the 
ODE mathematics becomes perceivably less transparent (cf., e.g., the classification 
problem solved for the first nontrivial double-singularity tobogganic model in [T6]). 
Indeed, under the presence of singularities in the potential V e ff(q) the wave functions 
<p n {q) generated by the corresponding Schrodinger Eq. ([T]) get analytically continued 
to a topologically nontrivial, punctured and multisheeted Riemann surface TZ. Even 
the description of this surface itself (characterized by its so called monodromy group 
[20] ) strongly depends on the details of the form of the potential. 

In the majority of the above-mentioned applications, the technical obstacle emerg- 
ing with nontrivial monodromies is most often circumvented via an ad hoc constraint 
which guarantees a "trivialization" of the topological structure of 1Z and of its mon- 
odromy [21J. An explicit form of the "trivialization" constraint is sampled, e.g., by 
Eq. (3) in Ref. [2] or by Eq. (6.31) in Ref. [7J. In this context, new an interesting 
questions emerge, of course, when one admits a nontrivial monodromy [T3J [T5] . 

The latter step of generalization is now in the center of our interest. In general, 
it would induce a rich menu of curves supporting the wavefunctions and connecting, 
in general, different sheets of the surface 1Z [IB]. In the spirit of an older letter 
[2"2"] . different spectrum may be attributed to each attachment of a new tobogganic 
curve to the same differential Schrodinger equation. In the two-singularity case, 
for illustration, the menu of these curves has been made available via a specific 
interactive demonstration using MATHEMATICA by Novotny [23J. 

In this context we feel curious what happens to the spectrum after the changes 
of topology in the most elementary though still nontrivial single-singularity case. 

3 Paths encircling single singularity 

For the sake of definiteness we may define our curves of integration, say, by formula 
Nr. (10) of Ref. [UJ, 

q(s) = q (N \s) = -i [i(s - ie)] 2N+1 , s G (-oo, oo) . (5) 



i 



Without any change of the spectrum these curves may be also smoothly deformed 
of course (cf. Figures 5 or 6 in Ref. [1] for illustration). In this case we only have to 
preserve both the left and right asymptotic parts of the curve q( N \s) unchanged (a 
detailed discussion of this point may be also found in [I]). 

From the point of view of physics (discussed in more detail, say, in [H1E32]) one 
of the most important features of the complexified quantum systems represented by 
Eqs. (CQ) + ([2]) on a line ([3]) (with e ^ 0) or on a curve ([5]) can be seen in the manifest 
loss of observability of their "coordinate" q R. In principle, the only measurable 
quantity remains to be the bound-state energy spectrum. This clarifies the main 
purpose of considering our specific class of a non-Hermitian models at nonvanishing 
winding numbers N = 1, 2, . . .: We may expect that the new free parameter N opens 
a new, potential-independent way of the modification of the spectrum needed, say, 
during its better fit to some given experimental data (of course, just hypothetical 
data in our present methodical study). 

In this broad framework the purpose of our present paper lies in finding a sample 
of quantitative analysis of the N— dependence of the energies. The two main sources 
of inspiration of such a study may be seen in the existence of a gap between the 
older, non-tobogganic results of Refs. [23] [reproduced also in Figure 8 of review 
[1] and using, incidentally, just £ = a = in Eq.(j2])] and the Bila's preliminary 
tobogganic numerical results which were presented in his very recent paper [19j and 
in his dissertation [25J. 

On the technical level the main ingredient of our present calculations will lie in 
the rectification transformation which will replace our imaginary cubic differential 
Schrodinger equation defined along one of the winding trajectories (JSj) by another, 
equivalent differential Schrodinger-type equation obtained via a suitable change of 
variables. 

In a more explicit specification of such a transformation we shall 

• distinguish between the non-winding curves where N = [abbreviating q^(s) = 
y(s) in Eq. and the genuine winding-paths q( N \s) with N = 1,2,... [use 
Eq. (jSj) and abbreviate q^ N \s) = z(s) at any N > 1 ], 

• start from the initial cubic and tobogganic bound-state problem 
d 2 1(1+1) 



dz 2 z 2 



^\z) = E^^\z), (6) 



• employ the change of variables z = — i (iy) 2N+1 or, in opposite direction, 
iy = (iz) a where a = 1/(2N + 1) is rational, i.e., where several (logarith- 
mic) Riemann sheets of the original variable z are mapped into the single (cut) 
complex plane of the new variable y. 

The detailed realization of such a project may be summarized as the direct use of 
the elementary substitution rules 

^-l^iV + l)^, -._J=tfL_- (7) 
<pW(z) = y N Mv) , L = (2N + 1) (i+D-l (8) 



2 J 2 



under which our original ordinary differential Eq. ([6]) becomes transformed into the 
following Sturm-Schrodinger [26] differential equation 

Mv) = ( 2N + !) 2 y 4N E P My) ■ 

(9) 

This change of variables realizes the one-to-one correspondence not only between the 
two respective (in general, logarithmic) Riemann surfaces of z and y but also between 
the original integration path (jSJ) (living on the original, tobogganic z— surface) and 
the new integration path for y which, by its present specification, coincides with the 
straight-line formula ([3]). 

During the search for the energies E[ N] the new, "rectified" , generalized eigenvalue 
problem (J9j) + (j3J) must still be addressed by purely numerical means in general. 
Still, the key merit of its choice is that one can employ some more or less standard 
methods. In particular, we intend to use here an appropriately adapted form of the 
large— i approximation method of Ref. [H], to be first briefly summarized in Sec. H] 
and subsequently applied to our problem in Sec. [5j 



dy y 



4 The large— i approximation technique 

The key ingredients of the approximations where the quantity l/£ remains small 
may be most easily illustrated via the non-tobogganic, single-sheeted version of our 
imaginary cubic oscillator, integrated along the straight-line path q 
Eq. ©, 



y(°) I 
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^ ] (q)=E^^(q). 



(10) 



The main idea of the approximative search for eigenvalues E n is based on the in- 
variance of the spectrum during any parallel shift of the path q(°'(s) pp. Thus, 
we are allowed to choose the free parameter e = Sj in Eq. (|2D in such a way that 
one and only one of the roots q = Q = Qj, j = 1, 2, . . . of the auxiliary equation 
dgV e ff(Q) = will lie on the corresponding specific, j— dependent exceptional line 
of our "coordinates" in Eq. ffTOl . q^(s) — > qj(s). 

In the next step we may rewrite our auxiliary equation in its fully explicit form 
2£(£ + 1) = 3iQ 5 and notice that in the vicinity of any of the quintuplet of its roots 
Q = Qji J = 1, 2, . . . , 5 we may Taylor-expand our effective potential in terms of the 
new, shifted complex variable £ = q — Q, 



1 



Veff(q) = V eff (Q) + - V: ff {q) £ 2 + - V% f {q) £ 3 + 



2 e ^ w/ ^ 6 

Next, we have to take into account the available explicit formulae for our V e ff and 
for all the roots 



Qj = —ir exp 



2ivr(j 



T 



1 



1/5 



J = l,2, 



'121 



Due to our assumption of the smallness of the quantities l/£ and/or 1/r we can 
choose j — 1 and simplify our Taylor series fTTTl) to its much more explicit form 



Veff(q) 



-- T 

2 



+ ^ 2 -5i£ 3 + 0(r- 



(13) 



f: 



This formula shows that up to the negligible error term our potential coincides with 
the deep and real harmonic-oscillator well, complemented just by a purely imaginary 
cubic perturbation. 

In the spirit of perturbation theory we shall accept our tentative choice of j = 1 
in what follows. The more general discussion of difficulties emerging in connection 
with the use of the other roots with j > 1 may be found in Refs. [H] and [27J. As 



long as we now have the new independent variable £ = (s) — Q\ = s which is 
real, we may insert the Taylor series (fT3j) in our Schrodinger Eq. (fTUl) and rewrite it 
in its leading-order anharmonic-oscillator reincarnation 



" J 5 3 + 15^2_ 5i ^ + o0 



d£ 2 2 2 



-iT + = 4V 0l (-ir + 0- (14) 



In the final step of our considerations we introduce another very small auxiliary 
quantity a = l/r 1//4 and rescale the coordinate £ — > cr£. This yields the final version 
of our approximate eigenvalue problem, 

" ^ + - 5a 5 if + 0(a 6 )] ^(-ir + = v 2 E^ ^](_ ir + . 

(15) 

As long as a <C 1 is very small, the leading-order version of this equation has the 
purely harmonic-oscillator form 

~ - + + ^)] ^ [0] (-- + = ^ Ei n ] ^ (-- + (16) 

which is exactly solvable. 

As long as we are interested here just in the very global features of our sample of 
tobogganic spectrum (like, e.g., the leading-order winding- number-dependence of the 
removal of its degeneracy, etc), we shall skip here the construction of higher-order 
corrections completely. Interested readers may find an instructive sample of such 
calculations in Ref. [28]. In this case we may already conclude that the low- lying 
spectrum evaluates to the following formula 



Ef =~ + ^ (2n + 1) + 0{r-^) , „ = 1, 2, . . . 

with an asymptotically vanishing error term and predicting the perceivably negative 
and approximately equidistant leading-order spectrum. 



5 The large— £ approximants at tobogganic N > 

On the basis of our preceding non-tobogganic = results we are now prepared 
to move to the genuine tobogganic models with windings N — 1,2, . . . entering our 
rectified Sturm-Schrodinger differential Eq. ([9]). Immediately, we may make use of 
the fact that one of the most convenient small parameters p = l/(£+ 1/2) 2 remains 
A^— independent. This allows us to keep the A^— dependence carried solely by the 
adapted shifts e = in the straight lines y(s) of Eq. (J3]). 

In the first step we have to Taylor-expand our effective potential again, 

V eff [y(x))] = ^_±!2 + i (_!)" {2N + 1)2 y ioN + s 



•7 



In contrast to section H] an amended notation must be introduced now since the 
original tobogganic (i.e., multisheeted) z— reference complex points Qj = gj; (0) 

differ from their transformed, Sturm-Schrodinger maps or partners Tj = q^°\o). The 
latter points lie on the rectified, straight-line contours ([3]) in the simpler complex 
y— plane with, possibly, just a single cut directed upwards. 

The positions of the latter reference points will again be determined by the 
derivative-vanishing condition d y V e ff(y) = 0. Thus, for our particular dynamical 
model each L— dependent root y = T = Tj(L) of this equation will coincide with one 
of the ION + 5 distinct roots of the following elementary algebraic equation 



2L(L + 1) = (2N + 1) 2 (1(W + 3) T 10N+5 . 

All of these roots form the vertices of a regular (10 AT + 5)— angle in complex plane. 
At N = we would just return to the pentagon which has been obtained in the non- 
tobogganic imaginary cubic model above. At any other iV > we shall abbreviate 



r = T ™ 



2L{L 



(2N + l) 2 (10N + 3) 



l/(10JV+5) 



(17) 



and set 

T x = -it , Tj = 2}_i e 2i7T/{10N+5) , j — 2,3, ... , ION + 5 . 

In a way which parallels our above discussion of the model with = we use just the 
same exceptional root T% = —it giving the maximal size of the shift e = r in Eq. (j3J). 
This choice also offers the single isolated real minimum of our effective potential along 
this line so that we may work again with the adapted real independent variable £ = s 
everywhere in what follows. 

After a somewhat tedious calculation one obtains the Taylor-series asymptotic 
approximation of V e ff(y) at any N, 



Veff (y) = -l(2N+l) 2 (10N + 5)T 10N+3 + cof N) t 10N+1 e + KN) T 10N ^ + O(t 10N - 1 ] 



with 

u (N) = (2N + 1) 



(ION + 3) (ION + 5) 
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Similar closed formulae can also be derived for fJj(N), etc. This means that our 
rectified Schrodinger equation will acquire the following approximate N— dependent 
asymptotic form 



^- 2 +u;f N) T^e+^N)T^e+ 



Vn(-iT + £)= (18) 



IeW (2N + l) 2 [t an + ...]+ 1 (2JV + l) 2 (10iV + 5) 



.1CWV+3 



</v(~ir + 



In a preliminary step of analysis we notice that the anharmonic-oscillator corrections 
proportional to the second coupling fi{N) will remain asymptotically negligible and 
that they may be ignored. Thus, in the asymptotic regime where r > 1 we shall 
apply the same considerations as tested in Appendix A. 



In an explicit error analysis we rescale £ — > cr£ where a = r ( 10iV + 1 )/ 4 . This leads 
to the equation 



d 2 



(19) 

which shows that although the asymptotic order of magnitude of the energies remains 
n— and iV— independent, En^ oc 0(£ 6 ^ 5 ), there is in fact no degeneracy in N and 
that also the degeneracy in n becomes removed by the next-order correction (cf. 
illustrative Figure [2]). 
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Figure 2: Tobogganic spectra fl20l) at p = l/(£ + l/2) 2 . The first five lowest levels are 

and displayed at the first four winding 
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In a more explicit asymptotic study of our model we have to keep in mind that 
our large auxiliary parameter r also changes with the winding number N, r = T(jv)- 
The source of this dependence does not lie only in the manifest presence of N in 
formula (jTTI) but also in the N— dependence of the upper-case parameter L. This 
is a mere abbreviation of the expression L = L(N,£) = (2N + 1)(£ + 1/2) — 1/2 in 
which one recognizes the presence of the small variable p = l/(£ + l/2) 2 used in 
Figures [2] or [31 Thus, the physics hidden behind the asymptotics in L can be re-read 
as a combined centrifugal-repulsion effect where both the coupling strength 1/p and 
the number of rotations iV appear to act in a coherent manner. 

At a fixed topological characteristic iV of the tobogganic quantum system in 
question our asymptotically simplified equations indicate that the level-splitting will 
predominantly be mediated by the harmonic-oscillator term £ 2 in Eq. (fl9|) . Thus, 
the leading-order equidistance of the levels (cf. Figure [3]) will only be broken, in the 
next-order approximation, by the cubic term 0(£ 3 ). An explicit estimate of the size 
of the latter correction is still feasible. From Eq. (FI9"|) one immediately deduces that 
the anharmonic-oscillator part of our potential will contribute by the corrections of 
the order O(r-( 107V+5 )/ 4 ). They may be perceived, say, as a higher-order modification 
of the coefficient uj(N)- 

We should add that another specific feature of the tobogganic problem emerges 



n 




Figure 3: The decrease of the level spacings E n N ^ 1 — En = G^ N \ 



in connection with the nontriviality of the right-hand-side "weight" term 

W{y) = {2N + lfy iN 

entering Eq. (jHJ) at ./V > 1. Fortunately, the scaling argument can be re-applied here 
to give the estimate of the first non-constant £— dependent correction to W(T) in the 
form 0(r~( 14iV+7 )/ 4 )). This correction is even smaller than the contribution of the 
anharmonic force so that it can either be neglected (on the present level of precision) 
or, whenever needed, calculated and incorporated, say, by the appropriately adapted 
systematic perturbation techniques as described in Ref. [Ti] . 

We may now return to Eq. ffl8l) and summarize our observations in the closed 
formula for the energies, 



10iV + 5 



r 6iV+3 + 



+ 



2n+l 
2N+1 



(10iV + 3)(10iV + 5) 



_7V+l/2 



0(t 



-(6iV+3)/4\ 



(20) 



giving the low-lying spectrum at n = 1,2, 



This is our final result. Once we 
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recollect that r oc £ 2 /( 10JV + 5 ) this formula reconfirms our above estimate that E n 
£ 6 / 5 , i.e., that in the dominant order the power-law exponent of I does not vary 
with N. In this setting the formula specifies the numerical factor and clarifies the 
mechanism of the removal of the degeneracy between topologically nonequivalent 
toboggans. 

In the subdominant, next-order level of precision the next numerical coefficient 
in the spectral formula will still retain an explicit form. The dominant part of the 
distance 

Q(N) = E [N]^ _ E [N) between levelg 

appears to have the power-law form 
with another constant exponent, G^ N ' oc I 1 ' 5 . This is well illustrated by Figured 
where we see how the differences between constants lim^oo G^/i 1 ^ 5 remove the 
topology-independence of the gaps, i.e., their leading-order degeneracy in N. 

This effect slightly weakens with the growth of winding number N > 1. It is worth 
adding that the function G^/i 1 ^ 5 of iV happens to have its maximum at N — 1/2 



i n 



so that the observed monotonicity in iV does not apply to the "anomalous", non- 
winding curve /£ 1//5 which, incidentally, lies somewhere in between /£ 1//5 and 

G^/e/\ 

6 Summary 

Our present paper offers the first quantitative confirmation of the hypothesis that 
the spectrum of quantum toboggans can be well controlled not only by the potential 
and, independently, by the suitably complexified asymptotic boundary conditions 
but, in addition, also by some suitable topological, non- asymptotic characteristics of 
the curves of the complexified "coordinates". These characteristics were shown to 
lead to the potentially useful changes in the spectra. 

For the sake of brevity we analyzed only the simplest nontrivial quantum to- 
boggans which are characterized by the presence of a single, isolated branch point 
singularity of ip{q) in the origin. In this case the non-equivalent toboggans can be 
classified by the mere integer winding number N. Of course, this does not mean 
that the analysis of the spectrum is easy. For illustration we may mention the re- 
cent diploma work [18] where, for several potentials, the stability and applicability 
of a few most common numerical methods has been shown restricted just to the 
non-winding dynamical regime with N — 0. 

More stable though still preliminary numerical results are available in Ref. [19] 
where the first encouraging samples of the real spectra of toboggans were presented 
using N < 2. Although this analysis remained restricted just to the smallest expo- 
nents M G (1, 3), an empirical confirmation of the reality of the energies at M = 2 
has been achieved. This seems to be in correspondence with Ref. [2H] where the 
exceptional role of integer exponents M in certain non-tobogganic models has been 
revealed. 

Intuitively one could expect that in the tobogganic cases the domain of a robust 
reality of the energies will move to the higher exponents M > M^ N >. Indeed, at 
iV = 0, the explicit numerical confirmation of the rightward shift of the boundary of 
the interval of the acceptable exponents can be found in Refs. [23] . Unfortunately, 
our present approach is unsuitable for an analysis of similar phenomena since, e.g., 
Eq. (j3J) controls the reality of the energies just via the insensitive sum of a given M 
with our very large \2i + 1|. 

On the side of the advantages of our present large— £ recipe we may mention the 
transparency of the determination of the unique and optimal line ([3]) of complex 
coordinates. In the leading order the use of this path reduces the dominant part 
of the effective interactions V e ff (q) to the solvable harmonic oscillator. In addition, 
whenever needed, the series of subdominant corrections can be evaluated in recurrent 
manner. In this sense the message of our present paper is non-numerical, showing 
that the low-lying levels of our toboggans can be given by closed formula and that 
they appear to vary sufficiently strongly with the variation of the topological winding 
number N. 
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Appendix A. The large— t approximation in an ex- 
actly solvable model 

Let us consider, in a test of the method, the exactly solvable example of the linear 
harmonic oscillator 

V^[q{s))] = ^^+u 2 {s-ief (21) 

in the strongly spiked dynamical regime, £ ^> 1. In a preparatory step of the analysis 
at large £ let us Taylor-expand this effective potential near its (complex) stationary 
point q = Q( HO > such that 

d g V e f f °\q) = . 

For our model every such a point is determined by the elementary algebraic equation 
£(£ + 1) = co> 2 T 4 with an easy solution, 



Q = Q, = (-iV' T , j = 1,2,3,4, r 



(22) 



All of these four candidates for the (in general, complex) point where V e ff(q) attains 
its (local) minimum also satisfy the sufficient condition since the real part of the 
second derivative of the effective potential remains Q— independent and positive at 
all q = Q since V? ff (Q) = 8u 2 . 

In the next step let us define the distances £ = q — Q from the minima and 
generate the Taylor series 

VefAq) = v eff (Q) + \v: ff {Q) e + \v:; f (Q) e+.... 

For our toy model this gives 

V eff (q) = 2oo 2 Q 2 + Auj 2 i 2 + i^ 2 e + 0{Q- 2 ) . 

Under the assumption that £ is large, i.e., \Q\ ^> 1, the insertion of this truncated se- 
ries would transform our bound-state problem ([1]) into its approximatively isospectral 
partner. 

The unrestricted freedom of the choice of the parameter e > in Eq. ([3]) enables 
us to set e = t and to let the line q^ (s) cross the point Qi = — i r. This is equivalent 
to the replacement of our initial Schrodinger equation (Til), at all the sufficiently large 
\£\ ^> 1 and/or r 3> 1, by the new equation 

~ - 2^ 2 t 2 + Alu 2 e + Oir- 1 ) tp n (-i T + = E n {uo) <p n (-ir + £) . (23) 

In the leading-order approximation this replacement leads to the new closed formula 

E n (cu) = -u v / (2£ + l) 2 - 1 + 2uj (2n + 1) + C(r 2/3 ) , n = 1, 2, . . . . 

It is rather amazing to see how well this £ ^> 1 estimate reproduces the exact low- 
lying spectrum 

E n {oj) = lo (An + 1 - 21) , n = 1,2, . . . ,n max , n max <£+l/2 
as derived, for our toy model ( f2~Tl) . in Ref. [30J. 



